Shear oscillations in the hadron-quark mixed phase 
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We calculate the torsional shear oscillations in the hadron-quark mixed phase of neutron stars 
whose structure depends strongly on the surface tension of the hadron-quark interface. It is shown 
that such frequencies become around ten times as large as those in the crust region, and those 
depend strongly on the surface tension. Additionally, we find that, with the fixed stellar mass, 
the frequencies of fundamental torsional oscillations in the hadron-quark mixed phase are almost 
proportional to the surface tension. So, with the help of the observation of stellar mass, one might 
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be able to obtain the value of surface tension via the observation of stellar oscillations. 
PACS numbers: 04.40.Dg, 97.60.Jd, 26.60.Kp 



I. INTRODUCTION 
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Neutron stars are formed in supernova explosions, which arise at the final moment of massive stars. The density 
inside the neutron stars can become higher than the standard nuclear density of pg ~ 0.17 fm -3 . Since such a high 
' density is almost impossible to be realized on the Earth, neutron stars might be the most suitable "laboratory" 
to see the properties of matter at high densities. One of the possible ways to see the properties of stellar matter 



could be the observations of gravitational waves emitted from the compact objects. Since the gravitational waves 
with high permeability will bring us the raw information of the wave sources, one can see the stellar properties via 
the observations of the gravitational waves. Then, the observations of gravitational waves enable us to collect the 
astronomical data, to reveal the basic properties of dense matter (e.g., [Jtl]), and to prove the gravitational theory in 
the strong gravitational region (e.g., p43)- Actually, the worldwide projects are going on to detect the gravitational 
waves associated with the astrophysical phenomena involving compact objects @. Another way to see the properties 
of stellar matter could be the direct observations of global oscillations of neutron stars. Via such observations, one 
would know the stellar mass, radius, and equation of state (EOS). This method is often referred to as asteroseismology, 
which is quite similar to helioseismology for the Sun. 

Unlike the gravitation waves, fortunately the observational evidences of the neutron stars oscillations have been 
detected, i.e., the quasi-periodic oscillations (QPOs) in giant flares emitted from the soft gamma repeaters (SGRs). 
Up to now, at least three giant flares have been detected, which are in the SGR 0526-66, the SGR 1900+14, and the 
SGR 1806-20. Furthermore, through the timing analysis of the X-ray afterglow in those giant flares, the specific QPO 
y—i . frequencies have been also found in the range from tens Hz up to a few kHz Since the central objects in SGRs are 
considered to be magnetars, which are strongly magnetized neutron stars, the discovered QPOs in giant flares could 
be due to the neutron star oscillations. In order to explain these QPO frequencies theoretically, a lot of numerical 
attempts have been done not only by the torsional oscillations in the crust of neutron stars but also by the magnetic 
oscillations [Tol - [20| . In addition, ascribing the observed QPOs to the torsional oscillations in the crust of neutron star, 
it became possible to reveal the properties of inhomogeneous nuclear matter in the crust [2l| - [24^ . 

The structure of neutron star is considered as follows: the ocean of liquid iron exists in the vicinity of stellar surface 
up to the density ~ 10 6 — 10 8 g cm -3 , subsequently the crust region exists from the bottom of the ocean up to the 
density of order of po, then the fluid core exists in the higher density region. In the most part of the crust, the 
nuclei could form a bec lattice due to the Coulomb interaction, while the existence of exotic nuclear shapes at the 
bottom of crust is also suggested in the recent studies [25T432I. According to such studies, with increasing density, the 
shape of nuclear matter changes from sphere (bee lattice |33| ) . to cylinder, slab, cylindrical hole, spherical bubble, 
and uniform matter (inner fluid core), which is collectively called "nuclear pasta." On the other hand, there are still 
many uncertainties in the core region. For example, the hyperons could appear under the beta equilibrium of nuclear 
matter when the density becomes higher than ~ 2 — 3po , or non-hadronic "quark" matter might exist in the innermost 
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stellar core [34j, |35j. Since the EOS could change due to the presence of these exotic components, the structure of 
neutron stars dramatically changes [H,|33]- Furthermore, as for the hadron-quark (HQ) phase transition, there are 
additional uncertainties such as the EOS of quark matter and/or the deconfinement mechanism. Anyhow the HQ 
mixed phase may emerge as a consequence of the Gibbs conditions, supposing that the HQ phase transition is of the 
first order 38]. The properties of the HQ mixed phase strongly depend on the electromagnetic interaction and the 
surface tension at the hadron-quark interface, which are called "the fine-size effects" [39|, Hoj. In fact, the appearance 
of the HQ mixed phase with "pasta" structures is subject to a balance between the surface tension and the Coulomb 
repulsion [37l l4l| , which is similar to the situation in the crust [30| and kaon condensed matter [42j . 

The one of the possibilities to distinguish the finite-size effects on the HQ mixed phase in the neutron stars is the 
observations of the astronomical phenomena. Actually, we suggested that possibility by using the direct observations 
of gravitational waves emitted [43| . Such attempts are very challenging and still the literatures are very few. In this 
article, we explore the torsional shear oscillations to see especially the dependence of the surface tension. For this 
purpose, preparing the neutron star models with the HQ mixed phase, we will examine the frequencies of torsional 
oscillations with the relativistic Cowling approximation. 

This article is organized as follows. In Sec. HI1 we describe the equilibrium of neutron stars and the adopted EOS. 
In Sec. IIII1 we show the equations governing the torsional shear oscillations and the boundary conditions to determine 
the eigenfrequencies. Additionally, the obtained spectra of such oscillations will be also shown. At the end, we make 
a conclusion in Sec. IIVI We adopt the unit of c = G = 1 in this article, where c and G denote the speed of light and 
the gravitational constant, respectively, and the metric signature is (— , +, +, +). 



II. NEUTRON STAR MODELS 



In this article, we focus on the nonrotating neutron stars with the HQ mixed phase. The equilibrium configura- 
tion of such relativistic objects becomes spherically symmetric solutions of the Tolman-Oppenheimer-Volkoff (TOV) 
equations. In this situation, the metric can be expressed as 



ds 2 = -e 2 * dt 2 + e 2A dr 2 + r 2 [dd 2 + sin 2 Odcj) 2 ) , (2.1) 

where $ and A are functions with respect to radial coordinate r. The mass function m(r) is associated with the 
metric function A as m(r) — r(l — e~ 2A )/2, which satisfies 

m' = 47rr 2 e, (2.2) 

where the prime on the variable denotes the partial derivative with respect to r, and e(r) is the energy density. The 
distributions of the pressure p(r) and metric function $(r) can be determined by solving the TOV equations; 

p' = -(e + P)&, (2.3) 

# = m + 4?rr3p . (2.4) 
r(r — 2m) 

In addition to these equations, one needs to prepare the EOS to close the coupled equations. 

According to [37} , we adopt the EOS including hyperons with the HQ mixed phase, which is taken into account the 
finite-size effects. That is, the EOS for hadron phase is adopted the non-relativistic Brueckner-Hartree-Fock (BHF) 
EOS, while the EOS for quark matter is assumed a generalized phenomenological MIT bag model, where we assume 
that u and d quarks are massless while s quarks have the mass of m s = 150 MeV as in [371 . |43| . Depending on the 
surface tension at the hadron-quark interface, the HQ mixed phase could become the non-uniform structure, where 
the structures of droplet, rod, slab, tube, and bubble are considered. The knowledge of the value of surface tension 
is very poor, but that value is theoretically estimated around a « 10 — 100 MeV fm -2 [44j, |45j. Since it is difficult to 
produce the HQ mixed phase with larger value of tr, we especially adopt a = 10, 20, and 40 MeV fm -2 in this article. 
At last, for the lower density region, the above EOS should be connected to the hadronic EOS proposed by Negele 
and Vautherin [46]. The EOSs adopted in this article can be shown in Fig. [T]in higher density region, where the EOS 
composed of only nucleons is also shown for comparison, and then the stellar properties can be obtained by solving 
the TOV equations with this EOS. 

Although there are not so many literatures about the shear modulus in the neutron stars, the formula about shear 
modulus in the neutron star crust has been derived in the zero temperature limit as 

M = 0.1194 x "*( Ze ) 2 ; (2.5) 
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FIG. 1: Relationship between the total energy density and the pressure. 

where m and Ze denote the ion number density and the ion charge, and a is the average ion spacing defined as 
a 3 = 3/(4.7rnj) [UEI!]. This formula is derived from Monte Carlo calculations with the assumptions that the shear 
modulus is averaged over all directions and the ion forms a perfect bcc lattice. In order to apply the expression of 
shear modulus (|2.5p for the droplet region in the HQ mixed phase, we consider that rij should be the number density 
of quark spherical droplet in the hadron sea, while Ze should be the total charge included in the quark spherical 
droplet. On the other hand, it is little known concerning the shear modulus in the other pasta structures, but the 
elasticity in such region is expected to be lower than that in the droplet region [49|, |50j. So, as a first step to see the 
behavior of torsional oscillations in the HQ mixed phase, we assume that /x = except for the droplet region in this 
article as in [H, [24|. Then, depending on the surface tension a, the shear modulus in the droplet region of the HQ 
mixed phase can be shown as Fig. [3J where the labels in figure denote the corresponding values of a in the unit of 
MeV fm~ 2 . 
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FIG. 2: Shear modulus /i in the droplet region of the HQ mixed phase as a function of energy density, where the labels denote 
the corresponding values of a in the unit of MeV mi -2 . 



III. TORSIONAL OSCILLATIONS 

The stellar oscillations in the spherically symmetric stars can be classified into two families with their parities. 
The oscillations with polar parity involve the density variation and stellar deformation, while those with axial parity 
are incompressible. Thus, one can examine the axial oscillations with high accuracy even if one neglects the metric 
perturbations by setting bg^ v = 0, which is known as the relativistic Cowling approximation. With this approximation, 
the restoring force of axial oscillations is shear stress characterized by shear modulus /i. Such oscillations are referred 
as torsional oscillations and can be described by a single perturbation variable y, which is corresponding to the 
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angular displacement in the 
fluid four-velocity as 



direction. The variable y is associated with the (^-component of the perturbation of 



5u* = e -*d t y{t,r)—-d e P e (cos 



(3.1) 



where Pi denotes the £-th order Legendre polynomial. The perturbation equation governing the torsional oscillations 
can be derived from the linearized equation of motion [51] 
can be written as 



Assuming y(t, r) = e luJt y(r), such a perturbation equation 
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e 2A y = 0. 



(3.2) 



Then, imposing the appropriate boundary conditions, the problem to solve is reduced to the eigenvalue problem. As 
in [22|, [2J, [5l|, we impose zero-traction conditions at the both boundaries of the HQ mixed phase. Here, it should 
be noticed that the torsional oscillations are also excited in the crust region, but such oscillations are completely 
decoupled with those in the core region. 

Before doing the numerical calculations, we can make a simple estimation about the frequencies of torsional oscilla- 
tions in the HQ mixed phase. Typical value of the shear modulus in the crust region is around /Lt ~ 10~ 10 — 1CP 9 km -2 
[22j . While, as shown in Fig. [2] since the shear modulus in the mixed phase becomes about 10 3 times larger than that 
in crust, the propagation time with the shear speed defined as v s = (/i/p) 1 / 2 becomes about 10 times smaller than 
that in crust. Thus, one can estimate that the frequencies of torsional oscillations in the mixed phase could become 
roughly 10 times as large as those in the crust region. 

Fig. [3] shows the frequencies of fundamental torsional oscillations with £ = 2 in the mixed phase as a function of the 
stellar mass for a = 10, 20, and 40 MeV fm -2 . Considering that such frequencies in the crust region are around tens 
Hz, as the above expectation, the frequencies of fundamental torsional oscillations in the mixed phase could be ten 
times larger than those in the crust region. Additionally, one can see from Fig. [3] that the frequencies of fundamental 
torsional oscillations depend strongly on the value of a. In practice, the frequencies for a = 20 and 40 MeV fm -2 are 
^40% and ^120 % larger than those for ex = 10 MeV fm -2 . Thus, if one would identify the observed frequencies as 
the torsional oscillations in the HQ mixed phase, one can probe the properties of such exotic structure. On the other 
hand, with fixed stellar mass, we plot the dependence of the frequencies of fundamental torsional oscillations on the 
surface tension in Fig. [4] From this figure, one can see the frequencies are almost proportional to a. Thus, with the 
help of the other observation of stellar mass, one could determine the value of a from the observation of the frequency 
of fundamental torsional oscillations. 
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FIG. 3: Frequencies of fundamental torsional oscillations with I — 2 as a function of the stellar mass, where the labels correspond 
to the values of a in the unit of MeV fm -2 . 

On the adopted stellar model, the torsional oscillations with different values of I also exist as well as the I — 2 
oscillation. In order to see the behavior of such oscillations, for the stellar model with a = 10 MeV fm -2 , we plot 
the fundamental frequencies with I = 2, 3, 4, and 5 as a function of stellar mass in Fig. [5J where the labels of otg are 
corresponding to the frequencies of £-th order oscillations. From this figure, one can observe that the dependence of 
frequencies on I is stronger than that on er, and that the dependence of frequencies on a might degenerate into that 
on I . That is, the identification of the properties of oscillation with using the observation of only one frequency might 
be difficult, even if one would know the stellar mass of the source. However, if the observation of several oscillation 
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FIG. 4: Frequencies of fundamental torsional oscillations with I = 2 as a function of surface tension with different values of 
stellar mass. 



frequencies would be successful from the source with the known stellar mass, to explain the observed evidence all 
together, one could be possible to make a constraint in a and to identify £, i.e., the degeneracy could be solved. 
Probably, with the development of observation technology, such observations will become possible and we will see the 
properties of the HQ mixed phase via the observation of the stellar oscillations. 
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FIG. 5: For a — 10 MeV fm 2 , frequencies of fundamental torsional oscillations with £ = 2, 3, 4, and 5 are shown as a function 
of the stellar mass. 

Finally, the frequencies of 1st overtone of torsional oscillations are shown in Fig. [6] Again, one can see that the 
frequencies in the HQ mixed phase becomes around ten times larger than those in the crust region due to the large 
shear speed in the HQ mixed phase. In addition to this, one can see the strong dependence of frequencies on the 
surface tension, where the frequencies with a — 20 and 40 MeV fm -2 become ~58% and ~85% larger than those with 
a = 10 MeV fm -2 . Compared with the fundamental oscillations, the dependence of frequencies of 1st overtone might 
not be so strong, but this difference could be still observable. Now, it should be emphasized that the dependence of 
frequencies of 1st overtone on a is different from that of fundamental oscillations, i.e., the frequencies of 1st overtone 
with the fixed stellar mass are not proportional to the surface tension (see Fig. [7]). So, via the both observations of 
frequencies of fundamental and overtone torsional oscillations, one could be able to make a severer constraint on the 
surface tension. 



IV. CONCLUSION 



In this article, we have considered the torsional shear oscillations in the HQ mixed phase, which can be composed 
of a non-uniform pasta structure depending on the finite-size effects. This type of oscillation is completely decoupled 
with the similar oscillations in the crust region. As a result, we have found that the frequencies of such oscillation 
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FIG. 6: Frequencies of 1st overtones of torsional oscillations with I = 2 as a function of the stellar mass, where the labels 
correspond to the values of a in the unit of MeV fm -2 . 
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FIG. 7: Frequencies of 1st overtones of torsional oscillations with I = 2 as a function of the surface tension with different values 
of stellar mass. 



in the mixed phase could be around ten times larger than that in the crust region and those frequencies depend 
strongly on the value of the surface tension for the hadron-quark interface. It is also found that the dependence 
of frequencies of fundamental oscillations is different from that of 1st overtone. Additionally, the frequencies of 
fundamental oscillations with the fixed stellar mass are almost proportional to the value of surface tension. Thus, 
with the help of the observation about the stellar mass, one might be able to determine the value of surface tension 
with using the observations of the frequencies of torsional oscillations in the HQ mixed phase. Furthermore, the 
resulting frequencies of fundamental oscillations in the HQ mixed phase are order of 100 Hz. This means that some 
of the QPO frequencies observed in giant flares, for example 150 Hz and even 626.5 Hz in SGR 1806-20 and 155 Hz 
m SGR 1900+14 @, might be associated with the torsional oscillations in the HQ mixed phase. In fact, the duration 
times for the oscillations with 150 Hz and 626.5 Hz observed in SGR 1806-20 arc quite long, which means there might 
exist a specific mechanism to excite such high frequencies for long time. The consideration of the torsional oscillations 
in the HQ mixed phase will be able to solve the puzzle for the theoretical explanation of the QPO frequencies observed 
in giant flares. 
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